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Abstract 

Electrons in shocks are efficiently energized due to the cross-shock potential, which develops because 
of differential deflection of electrons and ions by the magnetic field in the shock front. The electron ener- 
gization is necessarily accompanied by scattering and thermalization. The mechanism is efficient in both 
magnetized and non-magnetized relativistic electron-ion shocks. It is proposed that the synchrotron emis- 
sion from the heated electrons in a layer of strongly enhanced magnetic field is responsible for gamma ray 
burst afterglows. 
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I. INTRODUCTION 



Electron energization is usually considered as a secondary problem at heliospheric shocks, 
where most attention is paid to ion heating and reflection. In astrophysical shocks, however, these 
energized electrons emit the observed radiation, and are frequently the only source of information 
about the remote astrophysical process. Gamma ray burst (GRB) afterglow is believed to be syn- 
chrotron emission from electrons accelerated in the shock that develop during the interaction of the 
expanding ultra-relativistic plasma into the interstellar medium (ISM) l^. Estimates (e.g., Ref. Q 
and references therein) suggest that the required average energies of electrons reach a sizable part 
of the relativistic ion energy, and that the magnetic field in the emission region should be highly 
amplified, however, the origins of the electron heating and the magnetic field amplification remain 
poorly understood. In this paper we propose a single mechanism that accomplishes both, and is 
driven by the preferential deflection of electrons versus ions, when the former are lighter than the 
latter, by a local increase in the magnetic field. 

The mechanism of electron heating in heliospheric shocks is widely understood as follows yO: 
Electrons are decelerated more easily than ions, either by growing coherent magnetic fields in 
quasi-perpendicular shocks or by small scale magnetic structures in quasi-parallel shocks. The 
developing charge separation, however small it is, results in the build up of a cross-shock potential 
which is a substantial fraction of the incident ion energy. It is this cross-shock potential which 
decelerates ions when they become demagnetized in a thin transition layer of a quasi-perpendicular 
shock. In quasi-parallel shocks the parallel component of the magnetic field does not effect the 
ion motion along the shock normal, so that ions effectively become demagnetized just ahead of 
the transition. The same cross-shock potential which decelerates ions should accelerate electrons 
along the shock normal thus transferring energy from ions to electrons. The efficiency of the 
process is reduced by the electron drift in the magnetic fields, during which they lose energy by 
drifting down an electric potential. The final step of the process, electron thermalization, can be 
achieved by turbulent scattering following plasma instabilities. 

The mechanism of the prompt electron heating in steady state magnetized shocks is well-known 
electrons become demagnetized in the shock front if the ramp width is smaller than their 
convective gyroradius, or when the cross-shock electrostatic field becomes sufficiently inhomo- 
geneous to drag them across the magnetic field. In heliospheric shock these conditions are rarely 
satisfied since shocks are rarely this narrow. Moreover, only that part of the cross-shock potential 
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which cannot be eliminated by transformation into a de Hoffman-Teller frame [SD can be effec- 
tively used for electron energization. However, the profiles become steeper with the increase of 
the Mach number [16Q so that the conditions for demagnetization may be achieved more easily. The 
transition layer of quasi-perpendicular non-relativistic shocks consists of several distinct regions 
tM, the steepest magnetic field increase is a "ramp" (whose width is less than the ion inertial length 
k = c/ujpi, LUp^ = AimuC^ / nii) and a large magnetic overshoot (whose width is of the order of the 
downstream ion gyroradius). The overshoot height is found experimentally to increase with the 
increase of the Mach number [8]. The ratio of the ramp width to the ion convective gyroradius 
~ U cos 6 /{Vu/^ii) ~ 1 / M, where 6 is the angle between the shock normal and the upstream mag- 
netic field Vlu = eBu/irtiC is the upstream ion gyrofrequency, and M = 0.^1 ^pe is the Alfvenic 
Mach number. In perpendicular shocks the ramp width can be as small as l^. = c/upe [^. 



Theory of electron heating in quasi-parallel shocks has been developed s less elaborately, partly 
because of the lack of coherent structure in these shocks. Observations lioi imply that the domi- 
nant electron heating process is the same as in quasi-perpendicular shocks and appear to illustrate 
the importance of the DC effects of the coherent forces for the physics of electron heating in 
shocks. 

GRB-generated forward shocks in the ISM are ultra-relativistic T > 20. These shocks are 
parameterized by a = Bl/AnnumiC^'ju ^ 1 (this is written in the shock frame but is invariant). 
They are very high Mach number shocks, since the corresponding Mach number M = 1/a. 
Based on numerical simulations, it is widely believed that such shocks may be formed due to 
the development of Weibel instability jlll] into ion current filaments surrounded by regions of 
enhanced magnetic field. The filaments are elongated along the flow direction, with the magnetic 
field nearly perpendicular to the shock normal. The magnetic field around the filaments reaches 
nearly equipartition values but the magnetic filling factor is low. The width of a magnetic region is 
expected to be of up to tens of electron inertial length while the length of the region over which the 
surrounding magnetic field is high is determined by the ion scale. Although there is no gyration 
in these structures high magnetic fields at small scales make them play the role of a perpendicular 
magnetized shock front in what concerns electron energization. 

In this paper we suggest that differential momentum transfer to ions and electrons, typical for 
steady perpendicular shock and filamentary shock as well, results in the buildup of a strong poten- 
tial drop, comparable to the upstream ion energy. The electrons are demagnetized and receive a 
significant fraction of the original ion kinetic energy directly from the dc electric field. The accel- 
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erated electron energy is converted either into gyration energy (by the coherent magnetic field in 
magnetized shocks) or random motion energy (by small scale magnetic fields in Weibel mediated 
shocks) thus resulting in the coUisionless heating. In both cases a region of strongly enhanced 
magnetic field is developed in the shock front, where the heated electrons should efficiently emit 
synchrotron radiation. We show that although the details of the mechanism differ in magnetized 
and non-magnetized shocks, the underlying physics is very similar, and the eventual efficiency 
does not depend on the magnetization. GRB afterglows may be explained, at least in part, by 
radiation from these heated electrons. 

In proposing a mechanism for electron heating based on charge separation we do not mean to 
deny the existence of other mechanisms, e.g. decay and merging of magnetic islands, which can 
operate even with equal masses of both species. However, because the Weibel shock is otherwise 
required to "wait" for a bootstrap process in which electrons are heated by magnetic field but mag- 
netic field growth is limited by electron temperature |fl2||, we suggest that in the case of realistic 
mass ratios even a modest degree of charge separation can help to jump start the coUisionless 
shock process. 



n. MAGNETIZED SHOCKS. 



As will be seen below magnetized shocks are more restrictive in producing efficient electron 
heating, yet the basic features of the mechanism are typical for non-magnetized shocks as well 
(with suitable modifications). Therefore, we start our analysis with quasi-perpendicular magne- 
tized shocks. 

Relativistic shock propagating obliquely in ISM becomes nearly perpendicular in the shock 
frame, because of the Lorentz transformation, 9 shock = disAi/lu ^ 1 (here 7^ ^ 1 is the Lorentz- 
factor of the shock relative to ISM or, alternatively, the Lorentz-factor of the incident plasma flow 
in the shock frame). The de Hoffman-Teller frame, which has the velocity VutanO along the 
shock front, does not exist for Vu ^ c and tan^ > l/7„. In what follows we consider first a 
quasi- stationary perpendicular magnetized shock front where the fields are given by = B(x), 
Ex{x), and Ey = const. 
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A. Demagnetization conditions. 



The condition for the demagnetization by inhomogeneous is the statement that the acceler- 
ating electric field straightens the trajectory faster than the magnetic field bends it. The condition 
can be derived in the simplest way by approximating the inhomogeneous electric field with a lin- 
ear slope while ignoring the magnetic field variations in the electron equations of motion. Then 
the motion is described hy v — vq, x — Xq oc exp(At). Imaginary A (A^ < 0) corresponds to the 
particle gyration in the magnetic field (magnetic bending prevails) while and A^ > results in 
the exponential acceleration across the magnetic field, that is, demagnetization [4]. Relativistic 
generalization of the calculations in Ref.Q is straightforward (see Appendix |B]) and gives 

-l{l+l'vl/c')ie/m,)^>Ql, (1) 

where fi^ = eB /rrieC. If © is satisfied, electrons are efficiently accelerated across the magnetic 
field and acquire most of the cross-shock potential at the demagnetization region. The condition is 
local and cannot be satisfied in the whole shock transition layer, since —dE^/dx > is required. 
Thus, the electrons can be demagnetized while crossing a part of the magnetic inhomogeneity, after 
which they may return to be magnetized and the acquired energy is immediately converted into 
their gyration energy. Alternatively, electrons become demagnetized if the inhomogeneity scale of 
the magnetic field (1/ B){dB / dx) is smaller than the convective electron gyroradius c~^^l^^. 

The above demagnetization condition is derived in a simplified assumption that the magnetic 
field is constant. While this is not the case inside the shock, numerical analyses [|4|] have shown 
remarkable agreement with application of the non-relativistic version of ^ at the upstream edge 
of the ramp, and ([T]) should be considered an estimate. 

Demagnetization is required for an electron to utilize the cross-shock potential, otherwise elec- 
trons simply E X B drift, and the energy gain due to the potential (E^) is balanced by the energy 
loss because of the motion along Ey. Once the drift is substantially suppressed a net energy gain 
is achieved The energy gain is determined by the potential drop across the demagnetization 
region. When magnetization is restored no further energization occurs. The acquired energy is 
converted into the electron gyration energy where demagnetization disappears. Further collision- 
less "randomization" occurs through gyrophase mixing in the nonstationary and inhomogeneous 
fields of the shock front, thus resulting in the coUisionless heating [3]. Maxwellization is not 
required for the existence of the shock. 
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B. Magnetic structure and cross-shock electric field. 



For the purpose of description we consider a one-dimensional and steady shock. The basic 
equations of the two-fluid hydrodynamics for this shocks are given in Appendix lAl The cross- 
shock electric field can then be estimated using the momentum conservation: 

y^T^^ H = const, = {p^Vx)- (2) 

Here (...) means averaging over the distribution function and the summation is over both species. 
The discussion below is based on the basic picture justified by observations simulations [13], 
and theory u4\ . that the front steepening stops at the width much smaller than the convective ion 
gyroradius which ensures ion demagnetization inside the shock transition layer, and the assump- 
tion that this basic picture applies to relativistic magnetized shocks. As a consequence, ions are 
only slightly deflected within the transition layer (ramp) while almost all current necessary for the 
magnetic field increase is produced by electrons, which (partially) experience E x B drif. The 
latter allows one to estimate the electron velocity as Vy ~ {c/ Airnue) {dB ^ / dx) . Before electrons 
are substantially heated the magnetic force should be balanced by the electric force so that 

Ex ^ --^-^Bl enuA<f) ^ AB^Sn, (3) 
QTrriu dx 

where we have taken into account approximate quasineutrality and neglected the change of the ion 
density. For a <^ 1 even slight deceleration of ions causes strong enhancement of the magnetic 
field, which results in the development of the cross-shock potential which, in turn, further decel- 
erates ions. A spontaneous small enhancement of the upstream magnetic field causes exponential 
development of the magnetic field increase at the typical electron length scale (see below). The 
corresponding electric field given by ([3]). 

Upon crossing this narrow region of the magnetic field increase and potential development ions 
begin to gyrate. Assuming the gyrating ions to be a cold beam, it is easy to see that the momentum 
flow Txx in the particles is very small where the ions have gyrated by 90 degrees and are moving 
nearly perpendicular to the flow (x) direction. If the shock is to be quasi-stationary, this must 
be taken up by some combination of magnetic and electron pressure. For a weakly magnetized 
shock, magnetic pressure balance would imply a magnetic field far larger than that dictated by 
shock jump conditions. Electron pressure would require significant cross shock potential. The two 
quantities are connected by equation ([3]), so the argument implies both magnetic overshoot and a 
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large cross-shock potential. Since the shock may be unsteady, this argument does not constitute a 
rigorous proof of either, however, it shows that ion reflection is likely to cause extremely chaotic 
conditions in which pressure balance without strong cross shock potential and magnetic overshoot 
would seem to require implausibly fine tuning. 

In order to know whether electrons are indeed demagnetized one has to know the spatial profile 



of the shock. Two fluid hydrodynamics predicts 111411 that a perpendicular magnetosonic wave 
steepens down to the slope determined by the electron inertial length /g- Following the general 
principles of [Q], we seek nonlinear wave solutions that are asymptotically homogeneous, that 
is, n — ^ no, Vx — ^ vq, Bq, Vy 0, when x — oo. In this case Ey = vqBq/c. In the 

usual quasi-neutrality approximation charge-separation is weak throughout the wave profile 5n = 
{l/A'Ke){dEx/dx) -C n. Further derivation is given in Appendix O and results in the equation 

\ 1 d '-fe db 



uU NdxNdx 
_ (l + a)(b-l)-ab{b'-l)/2Pi 
l-a{b-l) 

where iV = n/rio = vo/vx,N = {l-a{b-l)){l-a{b^ -1)/2P^)-^ and 6 = B/Bq. The obtained 
expression is similar to those obtained previously for nonlinear stationary waves in pair plasmas 



[|15|1 . It is easy to see that the equations predicts the slope scale of /g = c^/Yej t^pe- The ratio lejre = 
{melvTLi)^!'^ ^pYeO'^l'^ <^ 1 for typical parameters of gamma-ray bursts. Therefore, electrons are 
expected to be demagnetized. It has to be understood, however, that the above small scale requires 
corresponding electron drift along the shock normal to ensure the current necessary to sustain the 
slope. Trajectories of demagnetized electrons are straightened along the shock normal and their 
drift is substantially suppressed, so that the ramp steepening does not proceed to scales much 
smaller than those required by the demagnetization condition. From the expression for and 
Eq. @) one can see that the amplitude of the magnetic compression reaches the values b ~ 
1 / ^/o- for strongly nonlinear structures in a low-cr plasma, in agreement with the estimates made 
independently earlier in this paper. 

To summarize, the basic points are the following: a) electrons become demagnetized if the 
typical inhomogeneity scale becomes smaller than the electron convective gyroradius Cy/^/cUpe 
or the cross-shock electric field slope is sufficiently steep to satisfy ([U), whichever happens first; 
b) the cross-shock electric field is related to the magnetic field as in ([3]), so that the potential 
increases with B^; c) the magnetic field, and hence the cross shock potential, increase to high 
values because magnetic pressure has to compensate the decrease of ion T^x as described by Q 



); d) large- amplitude magnetosonic waves steepen down to the scales c^/^/upe, as described 
by dH), which follows directly from the assumptions of electron drift and quasineutrality; f) the 
magnetic field in these structures increases io B / Bu ~ 1 / ^/o" before the singularity Vx = is 
reached; g) according to ([3]) the cross-shock potential is a substantial part of the incident ion energy, 
and h) the estimates above show that electrons have to be demagnetized (width is less than their 
convective gyroradius or ([T]) is satisfied). While not constituting a rigorous proof, these arguments 
show the plausibility and self-consistency of the proposed scenario of electron demagnetization 
by inhomogeneous cross-shock electric field and consequent heating. While the above scenario 
is described in terms of a monotonic magnetic field and potential increase across the ramp, it is 
likely that in real shocks the ramp itself breaks into sub-structures and the electron heating occurs 

n 

as a series of electric spikes [|16|] 



ni. NON-MAGNETIZED SHOCKS. 



Non-magnetized shocks are characterized by a very weak (or zero) upstream magnetic field 
so that the upstream convective gyroradii of both species exceed the system size and coherent 
magnetic braking is impossible. Weibel instability 111 ill produces magnetic filaments ahead of the 
main transiti on 1 17|] . Strong electron heating appears to be necessary for Weibel mediation at cr < 
rjiTe/niiC^Y [|l2| l where 77 is a dimensionless number less than unity. Otherwise, Weibel turbulence 
is predicted to be rather small scale and weak, so that ion scattering is relatively inefficient. Small 
scale magnetic filaments, where the magnetic field is aligned perpendicular to the shock normal in 
tubes or sheets, scatter forward going electrons more readily than ions (as does the perpendicular 
magnetic field enhancement for magnetized shocks), even if the electrons are not fully magnetized. 
Any deflection reduces the speed along the shock normal, thus the inflowing electrons are slowed 
relative to the inflowing ions by the turbulent field. The structure is no longer one-dimensional and 
stationary so that (IA5I) - (IA6I) are not applicable directly and (|AH) - (IA2|) should be used. However, 
assuming that electrons are scattered essentially randomly but small scale fields of the filaments, 
and neglecting ion scattering, one can average the equations over the perpendicular dimensions 
and time scales smaller than the ion transit time. Let us consider a single particle motion in the 
filamentary structure, taking the latter as given. The equations of motion read (for any species) 

= qEx + qx ■ {vtr X Btr), (5) 
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d 

—ptr = qEtr + qvx{x X Btr), (6) 
where tr denotes _L x. Here we assume that Etr and Btr are small scale rapidly (in space and 



time) fluctuating fields [11711 . while contains a global coherent electric field also. Denoting by 
averaging over rapid fluctuations, we assume that Etr = 0, Btr = 0, vtr = 0, but E^ ^ 0, 
Vx 7^ 0, E^tr 7^ 0, B'^tr 7^ 0, and v'^tr 7^ 0. In the lowest order approximation the particle flow is 
along X and scattering can be treated perturbatively: 

Ptr ~ {qEtr + qVx{x X Btr))T (7) 

where r is a characteristic "collision" time. Approximating ~ c, substituting ([7]) into ([5]) and 
averaging over rapid fluctuations, one has 

2^ 



d T 
Vx-rPx = qEx + —[x- {Etr X Btr) - B^tr] , (8) 
dx m7 

which is written for ions and electrons as well. Here we substituted (d/dt) — > v.j.{d/dx). 



Simulations lllin show that the generated magnetic field patterns are advected toward the shock 
front at speeds intermediate between the incoming plasma and the rest-frame plasma. In this case 
the electric fields are substantially weaker than the magnetic fields in the shock frame, so that the 



X ■ (Etr X Btr) term can be neglected relative to the last term which is nothing but the magnetic 
braking due to filaments. We now involve the smallness of r expected from the Weibel instability. 



The fastest growing modes have a scale length between the electron and ion inertial length 111211 . 
This means that in considering electron scattering, which we propose as a physical origin of charge 
separation, the ion scattering term which is oc r/mj is small relative to the electron scattering term 
which is oc r/mg. Thus, while the two terms in the right hand side of the electron equation ([8]) 
may be comparable for electrons, the last term is neglected for ions. Therefore, the ion motion is 
described by 

_ d _ — 

^i'^^Pi,x = (9) 

and for Vi^^ = c (negligible scattering of relativistic ions), one has 

cApi ^ = -eA(j), = - / Exdx, (10) 



in complete analogy with what happens to ions in a magnetized shock ramp: ions are deceler- 
ated by the potential which builds up due to charge separation caused by more efficient magnetic 
braking of electrons. 



With the same approximation, the electron energy changes as follows: 

^(meC^7j = -eE^Vx - eEtr ■ Vtr (11) 

qr —2 



~ -eE^v^ + e'-^[E^, - x ■ {Etr x Btr)]. (12) 

Unless the last term just happen to cancel the first term on the right hand side, the electrons acquire 
energy which is of the order of the potential drop e A0. Since this is the potential which decelerates 
ions, eA0 ~ mj7oC^, therefore, 

A(meC^7e) ^ eA0 ~ mi7oc^, (13) 

so that electrons acquire energy comparable to what the ions lose. Although we have not rigor- 
ously proved that this cancelation is impossible, we may note that in highly turbulent nonlinear 
environment the second term is likely to be a highly erratic function of space and time and it 
does not seem likely that its average would cancel the first term. That the first term should be of 
significant size is based on the fact that electrons are more easily scattered than the ions by the 
electromagnetic turbulence and this naturally results in the systematic charge separation during 
the early stages of a Weibel shock. The efficiency of energy transfer is higher than in magnetized 
shocks where only about a half of the potential can be acquired by electrons. This is because the 
electrons remain almost completely demagnetized throughout the whole region where ions decel- 
erate. Yet the electrons do not acquire all the momentum lost by ions, because of their scattering. 
Part of the momentum is transferred to the electromagnetic field. The pressure balance in this case 
takes the form 

^{PxVx) H ^ = const (14) 



Simulations llll . 



1711 show that filaments are convected by plasma and merge, so that both the local 
and average magnetic field density increase toward the shock transition layer. This is consistent 
with ([14)) : when approaching the transition the ion momentum decreases, as well as E'^tr (the latter 
because of the growth of the typical width of a filament), while B'^tr should increase. Similarly to 
what happens in magnetized shocks, magnetic braking of ions is necessary to convert the energy of 
the directed flow into thermal energy and decelerate the ion flow down to sub-relativistic velocity. 
As a result, the magnetic field is expected to achieve locally the equipartition values. This is also 
the region where the electron scattering by the magnetic field becomes strong. Once the electrons 
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and ions completely thermalize the magnetic pressure should drop to much lower magnitudes. A 
transient region of a drastic local enhancement of small scale magnetic field forms. 

Summarizing, all basic features found earlier in magnetized shocks (differential magnetic brak- 
ing, buildup of a potential and electron acceleration along the shock, magnetic field increase to 
equipartition values, conversion of the directed flow energy into thermal energy) are also present 
in non-magnetized shocks; in the latter, the latter local inhomogeneous magnetic fields play the 
role of the large scale magnetic background of the former. Respectively, the spatial scales of the 
corresponding "ramp" and "overshoot" are different and determined by ion gyroradius in magne- 
tized shocks, and by the filament merging in Weibel shocks. 

IV. SYNCHROTRON EMISSION. 

Having proposed that electrons acquire a substantial part of the incident ion energy due to the 
cross shock potential prior to entering a region of a strong magnetic field, we can now estimate 
synchrotron emission from this region. The main radiating region in magnetized shocks is the 
overshoot, behind which the magnetic field drops to low values. The radiating region in non- 
magnetized shocks should include the filamentary region before and behind the magnetic density 
peak as well. The below estimates are valid for magnetized and non-magnetized shocks as well. 
Let a shock propagate with the Lorentz factor 70 into interstellar medium with the density rij^m 
and magnetic field Bism, with a = Bf^^/Simism'm'iC^ ^ 1. In the shock frame the upstream 
density and magnetic field are n„ = Uismlo, = Bism- The electron energy in the overshoot 
is a fraction of the incident ion energy, that is, 7e = /i7o//^- The overshoot magnetic field is 
Sq/Stt = /2?^u?TijC^7o. The electron density in the overshoot follows the ion density which re- 
mains of the same order as the upstream density, rif, ~ 71^. At the lower end of the energy 
spectrum,the electrons emit synchrotron emission with the characteristic frequency and power (in 
the shock frame), respectively, ujm = {eBo/niec)'^l, Pm = {A/?,)aTC'^l{Bl/^n), where ut is 
the cross-section of Thomson scattering. In the observer's frame the characteristic frequency is 
^obs = lo^m, and the emission from unit perpendicular area becomes {dP/dS)obs = llPmNs, 
where A^^ = nj-o is the invariant surface density of electrons. Here r^, is the effective length of the 
radiating region. The observed frequency and emission per unit perpendicular area are 
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{dP/dS)obs = '2.(JTmiC^flf2^lnl^^lJ, ^Vo. (16) 

The largest uncertainty is in Tq since there is no satisfactory theory of the relativistic shock structure 
(neither magnetized nor non-magnetized). In a magnetized shock the effective overshoot width is 
determined by the ion gyroradius in the enhanced magnetic field, mjC^7o/ei?o, times number of 
ion loops necessary for gyrophase mixing. The maximum overall length is expected to be of the 
order of the ion downstream gyroradius or less, that is, r^, < f^miC^'jQ/ B^, where /s may be 
substantially smaller than unity. Respectively, {dP / dS)obs ~ lO^erg/cm^s ■ {nl/B3)'jlQf 1/2/3, 
where we normalized with the typical parameters for interstellar medium: rii = njs,„/lcm~'^, 
i?3 = Bism/^fiG. For a typical gamma-ray burst 70 = 10 — 30 several hours after the burst, and 
7io = To/ 10. In Weibel mediates shocks the overshoot width is determined by the ion inertial 



length llnh . In this case the enhanced magnetic field is strongly inhomogeneous, so that the ef- 
fective radiating width is = /a{c/ojS}^ where /2 and /^ together take into account the filling 



factor of about 10-15%. Simulations [|l7l] show that in Weibel mediated shocks the peak mag- 
netic density region is of the width of ~ 50(c/u;pi), but the region where B^/Svr ~ 0.1?7,„mjC^7o 
may be by an order of magnitude larger. The effective emission region may appear even substan- 



tially wider if the magnetic field decays as power law lllSh . Modestly estimating for these shocks 
/i ~ 72/4 ~ 10^, one finds {dP/dS)obs ~ lO^erg/cm^s ■ n^^^7jg. For the isotropic equivalent 
emitting area lO'^^'cm^ the total emitted power is P ~ 10^'^erg/s(n^/i?3)7[Q/f/|/3 in the mag- 
netized case and P ~ lO'^^erg/s ?7,^^^7Jq for non-magnetized shocks, emitted at the frequencies 
ujobs ~ lO^^s^ '^i^^Tio/i /2^^- hoth magnetized and non-magnetized shocks the magnetic field 
behind the overshoot drops down, B^ ~ Bo\fo. Respectively, the radiation frequency drops by 
the same factor, while the emission power drops by the factor 1/cr. 

This radiation from a thin region of enhanced magnetic field may be a significant fraction 
of the total afterglow emission. Consider the ratio of the afterglow from the magnetic region 
and from the entire downstream region. The fraction of the proper hydrodynamic time scale, 
Th ~ R/'yoc, that an electron spends in the effective overshoot region, is given by To/rh which 
is ~ /z'TujcTh ~ /sfTT-iC^lo/^BismR for the magnetized overshoot and ~ /4/ujpiTh ~ for non- 
magnetized shocks. The ratio of the magnetic energy density in the overshoot region to the average 
magnetic energy downstream is ~ 1/cr. Electron energies may remain comparable due to effective 
turbulent collisions. The relative afterglow outputs from the overshoot region and downstream is 
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then ~ f2To/Th(T, where R ~ cTobsll, Tots being the observer time, so that 



-^overshoot /2 /s ' 10 ^ (17) 



-Pdownstream (5/3 /xG) (Tob^/lO^ s)cr 

for a magnetized shock and 

-^overshoot /2 /4 ' 10 



> 1 (18) 



-Pdownstream nj^^ (To5s/10^ s)(T 

for a non-magnetized shocks. For realistically low a the emission power from the enhanced mag- 
netic field region formally exceeds the emission power in the rest of the downstream region. And 
the typical frequencies are much greater as well. 

The cooling energy 7c is given by the condition Pm{lc)ro/c ~ meC^'jc and therefore 
7c ~ 1^1 <JTllfonism.f2 which corrcspouds to the cooling frequency in the observer's frame 
ujc ~ 0Jobs{fJ^1c/ filoY ~ 10^'^s~^(_B|/f/ni/2/3) for the magnetized shocks, and much higher for 
demagnetized shocks, which means that radiative cooling does not affect the described processes. 



V. CONCLUSIONS 



We have shown above that efficient electron heating in relativistic coUisionless shocks can be 
generated by a cross-shock potential, developing because of the preferential deflection of electrons 
by the magnetic field. The cross-shock potential, which accelerates electrons across the shock 
front, is of the order of the incident ion energy, independently of whether the magnetic braking is 
caused by a coherent (for magnetized shocks) or small scale (for Weibel shocks) magnetic field. 
Deceleration of ions together with momentum conservation eventually lead to strong enhancement 
of the magnetic field in a small region of the shock front. This magnetic field enhancement ensures 
final thermalization of ions and electrons. Synchrotron emission from electrons from this enhanced 
magnetic field region seems to be able to explain the observed afterglow emission from gamma-ray 
bursts, within uncertainty of our knowledge of plasma parameters there. 

At scales below the ion gyroradius, the most likely scale for Weibel turbulence, differential 
scattering of ions and electrons by magnetic filaments can cause charge separation and strong 
electric fields in the shock plane as well as a long the shock normal, so it may be non-trivial to 
distinguish a systematic cross-shock potential from a purely stochastic electric field. Nevertheless, 
we suggest that a good way to test the idea of a systematic cross shock potential is to compare 
the electric field patterns for simulated pair shocks with shock simulations having a realistic ion 
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to electron mass ratio. The mechanism we suggest, which is based upon qualitatively different 
scattering of electrons and ions, works only for large mass ratios. For pair shocks, on the other 
hand, electrons and positrons can be separated by small scale magnetic fluctuations, but there is 
no systematic charge separation along the shock normal. If the systematic, cross shock potential 
drop for electron-ion shocks is comparable to the stochastic component, it would demonstrate the 
effect we are proposing. 
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APPENDIX A: TWO-FLUID HYDRODYNAMICS OF RELATIVISTIC SHOCKS 

Basic equations of two-fluid relativistic hydrodynamics read 



where s denotes the species, z, j, k = 1,2, 3, and Tjo and Tij axe the components of the energy- 
momentum tensor: 




(Al) 



(A2) 



Toi = {cpi). 



(A3) 



Tij = {viPj). 



(A4) 



Here ((...)) denotes averaging over the distribution function. 
In the one-dimensional stationary case the equations reduce to 



nvx = const, 
d 

-^Tix = nq{Ei + eijkVsjBk/c). 



(A5) 



(A6) 
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These equations should be completed with Maxwell equation with 

P = ^nsqs, (A7) 

s 

jk = ^ nsqsVs,k (A8) 

s 

APPENDIX B: ELECTRON DEMAGNETIZATION 

It is known [|4|] that in narrow nonrelativistic shocks electrons become demagnetized and effi- 
ciently heated due to the cross-shock potential. In order to find out whether such demagnetization 
is possible in relativistic shocks we reproduce the derivation of Ref. ^ with relativistic correc- 
tions. Namely, let us assume that a relativistic electron enters an inhomogeneous electric field 
Ex = {dE/dx)x, while the magnetic field inhomogeneity will be neglected. It has been shown i4\ 
that electron demagnetization occurs when two initially close trajectories diverge exponentially. 
Let us consider two close orbits xi{t), yi{t) and X2{t), y2{t), each of which is a solution of the 
equations of motion 

^{mv^'j) = -eE^ - eVyB/c, (Bl) 
d 

— (mwy7) = -eEy + ev^B/c. (B2) 

The equations for the differences 6x = X2 — xi, 6y = 1/2 — yi, 6vx = V2x — vix, Svy = V2y — viy 
can be easily obtained taking into account that ^7 = 'y^{vxSvx + Vy6vy): 

^[7(1 + -f\l/c'')5vx + -fHvxVy/c'')5vy] = - ^Svy, (B3) 

^[7(1 + -f^vl/c')6vy + i'^{vxVy/c')6vx] = n5vx, (B4) 

where we assumed for simplicity that B = const. Here = eB/mc and Ey = const. In the 
local approximation the obtained equations are linear equations with constant coefficients and 
substitution 5x, Svx, 5vy oc exp(At) gives 

[X'^^il + I'^vl/c') + ie/m){dEx/dx)]X-'6vx = -[fi + X^%VxVy/c^)]6vy, (B5) 
'jl/c^)Svy = [Q- X-f^{vxVy 

so that eventually 



A7(l + lh>Vc^)6vy = [n- X^'^{vxVy/c^)]6vx, (B6) 
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The local criterion of instability would read 



-(e/m)^>l]V7(l + 7W). 



dx 

For electrons entering the shock without gyration Vy 
divergence when 



(B8) 



e dEx 
m dx 



— >o, 
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and Vx ~ c, so that one gets trajectory 

(B9) 



with the divergence rate of 



A 



7 



-3/2 



e dE^ 



7 



1/2 



(BIO) 



The demagnetized electrons are accelerated by the electric field E^ across the magnetic field up to 
the point where the demagnetization condition ceases to be satisfied. At this point electrons begin 
to gyrated and all acquired energy is converted into their gyration energy. Beyond this point the 
only energy gain is due to the adiabatic conservation of the magnetic moment in the increasing 
magnetic field. 



APPENDIX C: NONLINEAR WAVES AND RELATIVISTIC SOLITON 

We consider a stationary perpendicular wave, d/dt = 0, d/dy = d/dz = 0, in the framework 
of the two-fluid hydrodynamics of cold relativistic electrons and protons (s = e,i for electrons 
and ions respectively): 



d 

' dx ' 
d 



^sVsx — {lsVsx) = qs{Ex + VsyBjc), (CI) 



msVsx-^ilsVsy) = QsiEy - VsxBjc), (C2) 

7. = (1 - vl/c' - vllc-r^l\ (C3) 

n^v^x = const, (C4) 

Ey = const, (C5) 

dB. 



dx 



-Atx ^ qsUsVsyl c = 4:7ie{neVey - niViy)/c, (C6) 



dEx V — \ 

— — = 47r qsUg = A-Ke^Ui - Ue). (C7) 

s 

It is worth mentioning that in the non-relativistic limit these equations have the solution in the 
form of the magnetosonic soliton lll4ll with the width ~ c/ujpf., where uj"^^ = ATme^ /rrie and the 
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amplitude depending on the Mach number. It should be noted also that ~ c/upe is the dispersion 
length of linear perpendicular magnetosonic waves. 



Following the general principles of 11411 . we are seeking for nonlinear wave solutions which 
are asymptotically homogeneous, that is, n uq, iij. vq, Bq, Vy 0, when x — oo. 

In this case Ey = vqBq/c. We shall consider weakly nonlinear waves in the sense that deviations 
from quasi-neutrality (charge- separation) are small throughout the wave profile 

5n=- < n. (C8) 

47re ax 

This assumption will be verified a posteriori. In this case = Ui = n ^ = Vex = v^, and 
nVx = riQVQ = const. Within this approximation we immediately get 

rrieJeVey + mi'JiViy = 0, (C9) 

B^ B^ 

nvoVxinie'je + fn^i) + ^ = nv^irrie + mi)7o + — ^, (CIO) 

nvo{me'ye + mai) + ^^^^ = nvo{me + mi)7o + (dl) 

(C12) 

where 70 = (1 — t'o/c^)"^''^, and further 



Vx 1 - - 1)/2I31 



{meje + rriiji) = {nie + mi)7o - ^ _ ° {h - 1), (C13) 

(C14) 



1-(t(6-1) ' 

where a = Bl/ATmQ{mi + me)c^ 70, (3q = vq/c, and h = B^/bo. It is easy to see that is a 
monotonically decreasing function of b in the range 1 < 6 < 1 + 1 /a. 
Using (|C9l ) and (|C6l ) one obtains 

(me7e + mi'ji)viy _ c dB^ ^^^^^ 
me7e 47me dx 

In what follows we shall make the assumption that the energy content in ions is always much 
higher than in electrons, that is, 

mi7i > me7e, (CI 6) 

so that approximately 

crrie'^edBz 

mi'JiViy = — — . (C17) 

Airen dx 
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Following the path outlined in the non-relativistic analysis [|l4h we substitute (IC17I) into (IC2|) for 
ions to obtain 

-'^x-ri-A — ^)-^ = 'VoBo-v^B^, (C18) 
ax V47re^n/ ax 



or, after normalization, 



NdxNdx 
(l + a)(b-l)-abib'-l)/2Pi 
l-(r{b-l) 



(C19) 



where = n/rio = vq/vx- 

In the asymptotically homogeneous point one has 



.2 \ J2 



where ,^ = 6 — 1^1. This point is unstable when [3^ > o-/(l + cr), in which case [|14|] the 
solution should be of a soliton type (non-periodic wave). The electron Lorentz-factor 7e cannot 
be represented as a function of b, so that (IC19I) cannot be converted to a quasi-potential equation . 
However, we can use the fact that 7e > to define a new coordinate dw = Ndx/'je, so that 

d\ ^ {l + a){b-l)-ab{b'-l)/2Pi 
^J^^^ = l-a{b-l) • ^^'^^ 

The derived equation is valid provided that the flow does not come to a halt, > 0, that is. 



b<bc= + 2(31 /a < 1 + 1/^^- (C22) 

When b increases the right hand side remains positive until b(b — 1) = 2/?q(1 + o-)/a. For 
cr ^ 1 and 7o ^ 1 this means that the sign changes when b = \f2jo ^ 1. At this point the 
denominator 1 — (t6 ^ 1. It is well-known \\% that there are no soliton solutions for cr ^ 1 in the 



pair plasma, where 7e = 7j = 1 — (t(6 — 1). For a soliton solution to exist 

(l + a)(&-l)-^K&^-l)/2/3g 
\, ^^ = ° ^^^^^ 

has to be satisfied for b^ < b^. Although complete analysis is impossible here it is likely that a 
soliton solution would not exist for too low a for the electron-ion plasma as well. 

For the analysis of the solution behavior it is sufficient to know that 70 < 7e ^ li{rni/me). It 

1 /2 

is easy to estimate the typical inhomogeneity scalce as 1^ ~ {c/ujpe)'^e ■ For cr ^ 1 (typical for 
gamma-ray bursts) the highest achievable magnetic field amplitude would grow as bmax ~ 1/ cr^/^, 
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thus ensuring strong magnetic compression. Since ab <^ 1 always, the electron current can be 
estimated as follows: 

UqCC 



where we assume that electrons remain relativistic: since becomes sub-relativistic, Vey ~ c. 
Then the typical length of the magnetic field variation is 

B 



(dB/dx) 



Bobil-ab/2) - aby2) (C25) 



where M = mi/nif,. The maximum length is achieved when b ~ and 1 — cr6^/2 ~ 1, where 

/ ~ c/ujpi. For smaller b ~ a/ ^/a, a ^ 1, the length / ~ a{c/ujpi), while for highest possible 
b ~ and 1 — cr6^/2 ~ and the length becomes / ~ {c/ujpe){MaY/'^ . 

It has to be understood, however, that the obtained expressions provide only an indication of the 
character of the wave steepening. Indeed, strong magnetic compression and narrow width ensure 
that ions behave nonadiabatically and begin to gyrate strongly in the vicinity of the magnetic field 
maximum. The ion gyration makes the cold hydrodynamical approximation invalid. Thus, the 
derived equation (IC19I) provides a satisfactory estimate of the wave profile only at the upstream 
edge of the shock ramp [|l4ll . which nevertheless is quite sufficient for physical conclusions to be 
made. 

Using (|C1I) one can find 



^^[("^hf + ^llD^l] = -e{mi-fi + ^ele)^- (C26) 



Taking into account the above approximation me7e <C nii'^i, and the expressions (IC13I ) and (IC14I ), 
one gets 



miC^-iQ Jf,^ \(3lJ V 1 - o-(fe - 1) 
For the above approximation the potential from the asymptotically homogeneous point to the point 
where d'^b/dw'^ = is easily evaluated as ecp ~ 0.5mjC^7o. 
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